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I. INTRODUCTION 



The interplay between contributions due to Feynman and hard-scattering scenarios |0] , the two basic mecha- 
nisms determining the large-Q^ behavior of hadronic form factors, is a crucial problem for studies of exclusive 
processes in quantum chromodynamics. To make a meaningful comparison of their magnitude, one should be able 
to calculate both contributions within the same self-consistent framework. To implement such a program, and 
especially in order to be able to take into account effects due to the nonperturbative part of the hadronic wave 
functions it is natural to rely on the light cone quantization which enables a Fock space representation for the 
hadronic current matrix elements. An important issue which has to be addressed when studying exclusive ampli- 
tudes in the region where both hard and soft QCD processes are important, is how to implement factorization, i.e., 
separate perturbative contributions from those intrinsic to the bound-state wave function |H,^,D . In the light cone 
quantization, the pion form factor can (in principle) be determined from 

where the summation extends over all quark/gluon Fock sectors which have a nonvanishing overlap with the pion, 
\l/„ are the corresponding wave functions, [dxidkij_]n is the relativistic measure within the n-particle sector and 
di = {1 — Xi) or —Xi depending on whether i refers to the struck quark or a spectator, respectively. Instead 
of having to deal with an infinite set of wave functions, which describe both the low and high momentum 
partons one can use the factorization procedure which allows to express high momentum tails of these wave 
functions in terms of suitably defined soft wave functions with nonvanishing support only at low momenta. In 
principle, this may be achieved by constructing an effective Hamiltonian which has nonvanishing matrix elements 
within the low momentum subspace. This effective Hamiltonian incorporates couplings between low and high 
momentum degrees of freedom present in the bare Hamiltonian, through a set of effective potentials connecting 
the low momentum states. The eigenstates of the effective Hamiltonian determine soft hadronic wave functions, 
\t'„. Because of asymptotic freedom, the removal of couplings to high momentum Fock states and construction 
of effective operators can be performed in a power series expansion in as- A similar procedure when applied to 
the current q'y^q yields an effective current operator. The contribution from the lowest order one-body current, 
as given by Eq. (|l|) is modified in higher orders by contributions from two-body and more complicated operator 
matrix elements resulting in 

F.,{Q^) = ^®T®^ (2) 

nm 

where summation over n, m now extends over the low momentum states only, and are the partonic matrix 
elements of the effective current operator. The dependence on the scale separating low and high momenta is 
indicated by /i. The matrix elements, T^n mix low momentum partons described by the wave functions, with 
intermediate high momentum states. The free energy of these intermediate states is bounded from below by and 
in order for the perturbative approach to make sense /x has to be chosen to be much larger than Aqcd so that 
tts = a{fj,) is small. At large momentum transfer we may ignore the dependence of Tmn on small transverse momenta 
kj^ and 1^ present in the soft wave functions. After that, (a;^ , k^ ; /i) and , 1^ ; /i) can be integrated over 

kj^ and 1^, respectively to produce the distribution amplitudes {e.g., for the pion) depending only on 

longitudinal momenta. To avoid large logarithms 0{as log((3^//i^)) emerging from the effective current matrix 
elements, we can set /i ^ |qj^| = \fQ^ , absorbing the logarithms into the evolution of the distribution amplitudes, 
4i(x, Q^). In particular, the latter enter the well-known formula for the pion form factor | 



{x,Q 
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■ dx 



(3) 



In the light cone gauge — 0, the distribution amplitudes are determined by matrix elements of composite 
operators in which the covariant derivatives coincide with the ordinary ones. In other words, to the leading 
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order in as and , only the minimal, valence Fock sector contributes. Conversely, any attempt to improve the 

leading order formula, for example by keeping the transverse momentum dependence of T„m, requires to take into 
account the presence of nonvalence degrees of freedom. 

Another important observation, supported by QCD sum rules |^,^ and by various constituent quark model 
studies , is that the soft contribution to the form factor may be quite substantial. To make a self-consistent 

comparison of the relative size of the purely soft term and those containing the hard gluon exchange, one should be 
able to calculate them using the same formalism, which is not necessarily based on the expansion but which 

rather enables to compute the perturbative corrections to the soft form factor. In other words, we propose to take 
advantage of the formula given in Eq. (|^), retaining there the soft wave functions, rather than expressing the form 
factor in terms of the distribution amplitudes as given in Eq. (|^) for the leading contribution. 

In practice, to develop a 3-dimensional renormalization program for a Hamiltonian which also produces light-cone 
bound states is rather difficult. One of the reasons is that cut-offs violate Lorentz and gauge symmetry, and thus 
it is highly nontrivial to achieve the necessary cancellation of IR divergences. 

It is also unclear how one could systematically incorporate the pQCD radiative corrections into a bound-state 
approach with phcnomcnological low-energy interaction kernels. 

A simple approach which enables the connection between the three-dimensional soft wave functions and the high 
energy scattering kernels of perturbative QCD was proposed in Ref . jl^] where it was also applied to describe the soft 
and hard contributions to the heavy-to-light meson transition form factor. The method relies on Green's functions 
just like the QCD sum rule approach Let us recall that the QCD sum rule method uses the quark- hadron 
duality relation 

J^^ ph-<i-(.)e-/^'^^ ds = p<^"-k(.)e-/^^^ ds + Y: (j^, (4) 

between the hadronic spectral density p^^'*''°"(s) on one side and the perturbative spectral density p'^'^^''^{s) and 
the condensates (O^) on the other. The exponential weight e^"!^^ in this relation results from the Borel trans- 
formation which both emphasizes the lowest resonance contribution into the dispersion integral and improves the 
convergence properties of the condensate series generated by the operator product expansion. For an appropriate 
choice of the Borel parameter M, the left-hand side of the duality relation is rather closely approximated by the 
lowest resonance contribution while the right-hand side is essentially given by the perturbative term alone. In such 
a situation, one can approximate the lowest resonance contribution by the Borel transform of the relevant Green's 
function. Such an approximation may be rather accurate when the Borel parameter is taken within a specific range 
for which the higher state contribution into the left-hand side is close in magnitude to the condensate contribution 
to the right-hand side. As we will see below, such an ansatz is analogous to the exponential (oscillator-type) model 
for the soft pion wave function. Calculating the as(Af^) expansion for the Green function (i.e., higher-order terms 
in p'i"'^'''^(s)), we can construct a unified self-consistent model which includes both the lowest Fock state contribu- 
tion and the higher Fock terms generated by the radiative corrections. The model, e.g., automatically preserves 
such an important property as gauge invariance. The latter is crucial for cancellation of infrared divergences. 

In application to form factors, as described in Ref. ||l3|| , the Borel transformation is used as a method which, 
on a loop-by-loop basis, enables to explicitly identify the wavc-function-like contributions to the form factor which 
have the structure of the Fock-space decomposition Eq. (j^) . An important advantage of this method over a naive 
summation of perturbative corrections to model wave functions, is that no ad hoc assumptions about how to connect 
soft and hard regimes have to be made. The formalism constrains how the IR divergences are distributed among 
the wave functions and scattering kernels to guarantee the IR finiteness and factorization scale independence of a 
matrix element. As we will see, this requires that both the higher Fock space contributions, with, e.g., n — qqg^ 
and the radiative corrections to the current matrix elements, T„m should be present. 

The paper is organized as follows. In the next section we discuss the model and details of the form factor calcu- 
lation. The Sudakov suppression of the electromagnetic vertex is discussed in Sec. III. Conclusions are summarized 
in Sec.IV. 
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II. MODEL SPECIFICATION 



As an illustration of the method, consider a two-point amplitude defined by 

ip+fllip') = -i I d4:re^^'-(0|r{g(x)7+75(z(x),g(0)7+75(?(0)}|0). 



After the Borel transformation, 

n(/) ^ n(/3) = ^ jdp'e-P"'^P"li{p') = j dse-'/^^"p{s), 
where p{s) = (n(s — ie) — n(s + ie))/27ri, the bare quark loop contribution to 11 yields. 



jjquark^^^ = f d8e-'/^'^^p'i'''''^{s) =4x2x3 / [dXidl!ii_L]2 Cxp --L 

J J l •'^i 



(5) 



(6) 



(7) 



with the numerical factors in front of the integral coming form the trace over Dirac indices, the quark current 
matrix element u(A)7+75u(A') = 25\-\i, and sum over three colors respectively; the measure is given by 



The quark spectral function, pq^^'k jg gjyg^ ]-,y 

pquark(^-) = 2x/6y [dxd^kj_]2 2^/6 6 (^s - ^^^^^ ■ 
The contribution of the lowest hadronic state to n(/3) is given by 



(8) 



Ml 
2/32 



(9) 



(10) 



where in our case = is the mass of the ( = ) ground state, the pion, and fh = f-n is the pion decay 
constant 



(0|g(0)7+75?(0|P) = ihP+ . 



Therefore, if we assign 



T / 1 ^ 2V6 
W2(a;i, fcjj_) = -— cxp 



1 

2^92 



then the decay constant /^r may be expressed in terms of this valence wave function: 

= 2^6 j [dxid\^iA_]2'^2{xi,kiA_). 



(11) 



(12) 



(13) 



This expression agrees with the one obtained using the light cone quantization directly. In the following we will 
drop the subscript on the wave function and on the light cone measure when referring to the valence qq component 
of the gromid state wave function. In the valence sector, the individual quark momenta Xi,\Li±_ may be replaced 
by the Jacobi variables defined by a; = a;i, 2:2 = 1 — a;i and kj_ = kij_ = — k2_L such that 



[dxidk.i±_\2 — [dxdkj 



rf.xd2k_L 
167r3 



(14) 
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In addition to the Borel transformation we will use also another smearing procedure corresponding to the "local 
duality" prescription. Its origin can be briefly explained in the following way. In QCD sum rule calculations one 
typically uses the following ansatz for the hadronic spectral function, p''^'*''°"(s) 

phadro„(^) « flS{s - Ml) + p'i"^^\s)9is - so), (15) 

with the effective continuum threshold sq fitted from the relevant QCD sum rule Eq. by requiring the most 
stable [M^-independent] result for physical quantities. In many cases, the stability persists in the whole large-M^ 
region. Fitting sq and then taking the limit Af — > oo for the Borel parameter in Eq. (^) gives [ p^ 

r p^--r^i^s)ds. (16) 

JO 

The local duality states that the two densities, phadron ^quark gj^g ^y^q same result provided the duality interval 
So is properly chosen. Comparing Eqs. (|^) and ( |l^ ) results in 

f,^2V6j [dx^d^k^±]2■9^^ix^,k,^), (17) 

where the "local duality" pion wave function ^''^^ is defined by 

The parameter sq = 47r^/^ ~ (830MeV)^ serves here as a transverse momentum cutoff and plays a role similar 
to the wave function scale P in the approach based directly on the Borel transformation. In the following, we 
will continue using the Borel transformation when deriving analytical expressions. But we will also give numerical 
results obtained for the local duality wave function. 

A. Soft contribution 

To obtain a model expression for the soft contribution to F^^iCp) in a form of an integral over the valence light 
cone wave function, consider the three-point function 

{p+)^T{p\p'^) = - [ d4xd4ye^''--'f'^(0|rg(x)7+75g(x),g-(0)7+g(0),g-(2/)7+75g(y)|0) (19) 



in the kinematical region defined by — p+ = (7+ = and {p' — p)^ = = = —Q^- The quark loop 

contribution to the double Borel transform ation, 

T(p2,p'2) ->T(/3) - [ dsds'e-^'^^'e-'''^^" p{s: s), (20) 



IS 
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p^--\s, s') ^ SpT{s, s') = 2(2^6)^5 j^s - ^ f ) ^ ^ I^li±i^Sii_^ j , (21) 

and therefore the qq contribution yields 

ri-*(/3) = 2/1 exp J [dxdk^]^{x,k^)^ix,k^ - xq^). (22) 



The ground state hadronic contribution to T is given by 
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T''~{P) = 2f^F^iQ')exp 



2M^ 

2/32 



(23) 



Identifying Eq. (|2^) and Eq. 

to F^iQ^), 



3) yields the standard expression for the valence quark wave function contribution 



(24) 



As illustrated on the two examples above, for a single quark loop we have generated expressions for hadronic 
matrix elements in terms of model light cone wave functions. As will be shown below, this remains to be the case 
beyond the single quark loop. When the higher loops are considered, the nonvalence wave functions emerge. 

In the following we shall work in the chiral limit and set = 0. With the local duality wave function given 
by Eq. ( p^), the normalization condition for the form factor, F^(0) = 1 is satisfied automatically. Taking the wave 
function (|l2|) based on the Borel transformation, produces only a half of the necessary form factor normalization 
at = 0. In the standard light-cone approach, the remainder comes from contributions due to the higher Fock 
components like qG . . . Gq corresponding to the "intrinsic" gluon content of the hadron. The intrinsic gluons should 
be contrasted with those emerging from explicit radiative corrections to the two-body wave function. Since the 
local duality wave function provides 100% of the = normalization, one can interpret it as an effective two-body 
wave function absorbing in itself the soft part of the higher-Fock states. In the spirit of effective wave functions, 
we will also modify the Gaussian wave function and instead of Eq. (|l^) wc will use 



Mi 



1=1 



with the extra normalization constant N and the parameter /3 fixed by Eq. (Jl 
iV = 2, /3 = tt/^ - 400 MeV. 



(25) 



and -P7r(0) — 1 which lead to 



B. One-gluon exchange and hard contribution 

The one-gluon-exchange diagrams are shown in Fig. Ib-c. In addition, to order Us there are three self-energy 
type diagrams from dressing of the propagators of the triangle diagram shown in Fig. la. The hard-gluon-exchange 
contribution to i^^ will be defined as an IR-regular part of two diagrams (Figs. lc,d) corresponding to vertex 
corrections to the axial vector currents. The IR finite part of the diagram with electromagnetic vertex correction 
(Fig. lb) contains terms which produce the Sudakov suppression after all-order summation. Since the operators 
which define the correlator T correspond to conserved currents, the sum of all six amplitudes is UV finite. It is 
also IR finite since there is no net color flowing into or out of the triangle. 

Just like in case of the one-loop triangle diagram we calculate Feynman amplitudes corresponding to each two- 
loop diagram by performing analytically the integrals over the "minus" light cone components of loop momenta. 
Then, we apply the Borel transformation, and the net result may be schematically written as 

The first term comes from the one-loop diagram. The remaining four terms represent the 0{as) contributions. The 
T-term corresponds to the vertex correction to the electromagnetic current convoluted with two qq, i.e., valence 
meson wave functions. Two other amplitudes, Tgi and Tg2, are convoluted with the valence wave function on one 
end and the nonvalence qqg wave function on the other one. Finally, the J^g-term corresponds to the nonvalence 
qqg meson wave functions induced both in the initial and final states. 

As an example consider the diagram shown in Fig. lb. The Feynman amplitude is given by 



xTr 



(27r)4 (27r)4 

7+75 (# - «7"(/ - V)l+{i> - V)l^{i> ~ «7+75(-« 



{p' - kY{p' - iY{p - lyij) - fc)2fc2(fc - z)2 



(27) 
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The k and I integrals pick up poles in the spectator quark and gluon (l/(fc — l)^) propagators leading to 

2 f dxdk± dydl± y[l± - q±] ■ [l± + (1 - y)q±] 



1 



167r3 167r3 - yqj^]^ 

1 



1= 



i{l-x) + P x{l-x) y{l-y'){l-x] + 

1 1 



P' 



12 _ (kj_-a:qj_)^ 
x{l — x) 



+ ie p 



12 _ (kj_-a;q^)^ _ (U-i;q^)^ 
x(l-x) 



(28) 



After the Borel transformation, the contribution to i'V from this diagram can be expressed in a form of Eq. (|2£ 
Explicitly, 



dxdk I dvdl \ r , , . , , , , . ^ „p , , 

, , „ 3 i±)Tgi-^ix, - xq^) + ^{x, k^)T^2'^g{x, - xq^; y, 1^ - yq±) 

+-^g{x,kjL;y:l±)Igg'^g{x,kjL - a:^q±;y,U - Z/q_L)] , 



where 



dydl^ 



Niy,h 



167r3 y{l - y) D{y, l±)D{y, lj_ - yq±) 



(29) 



(30) 



and we have defined 



ySTTCpas N{y, L) p 
yil-y) D{y^^)' ' 



92 



V8^6Va, iV(2/,U) 

y(i-y) £*(?/, i±-yqj 



iV(2/,U) 

2/(1 - y) 



y(i - y) 



y(i - y)' 



(31) 



(32) 



The valence wave function 5* is given by Eq. ( psf ) (with A/^r = 0). The contribution from the qqg state in the 
electromagnetic vertex loop is determined by g which is given by 



*.(x,k,;,,l,) = ^^^— exp 



1 

2^02 



x(\-x) y(\-y){\-x) 



(33) 



Here y represents the fraction of longitudinal momentum of a quark carried away by the gluon. The spectator 
quark has longitudinal momentum x. The transverse momentum 1^ is the relative momentum between the struck 
quark and the gluon. Note, that just like the argument of the valence wave function ^'2 is determined by the 
invariant mass of the free qq system, the argument of ^'g is given by the invariant mass of the free qqg system. The 
contribution proportional to the overlap of valence wave functions is proportional to the vertex form factor, . 
The UV divergence of is canceled by a half of the sum of self-energy corrections to two quark lines connected 
to this vertex. The total contribution to from these two self-energy diagrams can also be cast in the form of 
Eq. ( p6| ) with Igg = and is given by 

F^{Q')^F^m= I ^^^{x,k^)T^^{x,k^-xq^) 

dxdk J. dydl± ^^^^^^ l^)T^^^^>{x, k± - xq±) + ^{x, k±)T^2'^g{x, k± - xq±; l± - yq±)] , 



(34) 
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with 



and 



IGtt^ 2/(1 - y) 



y 



D{y,l^) D{y,l^-yq^) 



y 



2/(1 -zj) DiyM) 



Tf2 



y 



2/(1 - y) D{y,l± - yq±) 



(35) 



(36) 



The remaining UV divergence coming from these two self-energy diagrams cancel the UV divergences resulting 
from diagrams which dress the axial current vertices. The IR divergence of corresponding to either |1_l| ^ or 
|1_L — yc[±\ — » is removed by two contributions in Eq. ( |29| ) proportional to the qqg wave functions. When both 
denominators in vanish, i.e., \l± \ and — yq± \ — > 0, all particles in the loop dressing the electromagnetic 
vertex go on-shell and this singularity is removed by the contribution proportional to the qqg wave functions in 
both initial and final state. Similarly, the IR divergent contribution to in Eq. ( ^5|) corresponding to either 
—^0 or |1^ — 2/q_L| — > is canceled by mixing with the qqg wave functions. 
Consider now the third and fourth diagram shown in Fig. 1. After integrating over the "minus" components of 
the loop momenta they can also be written in the form given by Eq. (^6|). In particular, in the limit |qj^| » f3 the 
dominant contribution comes from the first term in Eq. (p6[), i.e., the overlap of the valence wave functions and it 
is given by 



f.(q2)^f^ _;^(q2)^ j ^*^.i,(2^,i^)T_h(y,U;x,k^)^(x,k^-xq^), 



(37) 



where 



-^exch 



y-x 



x(l - x)\l + 2/(1 - 2;)ki + {yx + (1 - y){l - x))\^ ■ 



(l-L-k_L)^ 



1^ 



y{l-y) x{l-x) 



(38) 



To leading order in 1 /q^ this would naively give 



l&TTCpas f dxdk± dydl± ^{y, l_L)*(a:, k^) 



xy 



and therefore the expression for -FV.oxch becomes identical to the leading order formula, Eq. 
amplitudes given by the asymptotic formula 



c^{x) 



^*(,,k,) = ^/.x(l-x) 



(39) 

with the distribution 
(40) 



The second equality above follows from the normalization condition given in Eq. (^^ . Of course Eq. (|39|) does not 
exactly follow from Eq. (|3^) since the latter, and therefore also i^7r,cxch, is not well defined; Texch has IR divergences 
coming from two energy denominators in the square brackets in Eq. (pq). These singularities correspond to open 
qq and qqg intermediate states. These singularities are canceled in Eq. (Bq) by contributions from 5* ® Ti,2g ® 
and (g) Igg (g) r csp cct Ivcly. From Texch we may however define a sc/ieme-dependent hard contribution kernel, 
-^exch (/^) ^'^d the corresponding hard form factor from one gluon exchange. 



phard 
TT, cxch 



dxdk± dydl± 



*(2/, li_)T^^i^{y, 1±; x, k^; fj.)-i'{x, - xq^). 



(41) 



by cutting-off light cone energy denominators such that Eq. ( pT\ ) iVc5«:h becomes IR finite. For example, we may 
define 



^cxch (m) 



ne(A 



■exch 



rrihard 
cxch 



ne(A 



cxch; 



(42) 
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where Di is either one of the two denominators in Eq. ( p8| ) and Q{D) ~ 1, in the momentum region in which 
\D\ > /i^ and Q{D) for \D\ << /x^. The hard contribution from one gluon exchange, T^^^^f, defined in this way 
is IR finite but cut-off, /i dependent. The /i-dependent IR singular piece of Toxch, T^^y^ifJ') when combined with the 
contributions proportional to the qqg wave function produces IR finite but low momentum-dominated contribution. 
Of course, when this term is combined with the contribution from Tcxch (a*)' ^^"^ /i-dependence disappears. If we 
choose the factorization scale fi to be comparable to that setting the width of the soft wave function, i.e. fi ^ f3 (or 
~ So for the local duality wave functions) then contribution from nonvalence sectors becomes reduced by the 
residual contribution from T^f^j^ = Texch — T^c'^Sf , i-e. the IR dominated piece of the valence contribution. Then 
for |q_L| >> /i bulk of the form factor will come from the hard gluon exchange given by T'e'xch ■ With Q given by 



\D\ 



(43) 



the effect of the nonvalence contributions is to effectively add a mass term of order /i to the two energy denominators, 
Di in Eq. (pi) which leads to a similar effect as the gluon mass regulator used in Ref. |0 . 



III. SUDAKOV SUPPRESSION 

Before proceeding to the analysis of the numerical results, let us discuss how the Sudakov suppression from the 
gluon radiation of the struck quark results in our formalisms. The exchange of the gluon across the electromagnetic 
vertex, Eq. (pq) gives for the quark spectral function, p'^""'^(s, s') 



Pvert[S, S 



dxdk^ dydl± y[l± - q±] ■ [l± + (1 - y)q±] 



6 \ s 
S I s' 



x{l — x) 
{kj_ ~ x(\x_\ 
x{\ — x) 



IGtt-^ Voir" 

-s(s- 



1^ 



x{l-x) y{l~y){l-x)^ 

(k_L - xqi)2 (U - 2/q^)2 



x{l-x) y{l-y){l-x) 



(44) 



We now write the first term in brackets as 



S s- 



x{l — x) 



r:(l-x) 



11 



11 



k2 



Xil -x) Ij , 

v(l-V)(l-x) ~^ x(l-x) 
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1^ 



x{l-x) y{l-y){l-x) 



(45) 



and similarly for the second term with s —> s' and kj^ — > kj^ — a;qj_, 1_l — > 1_l — ?/q±- In the first term on the 
right-hand side in Eq. (Bsh the qqg invariant mass, \\/{y{\ — j/) (1 — x)) in now regularized in the IR by the incoming 



qq virtuality, \c\_/x{l — x). Thus when integrated in Eq. (|4J) over gluon transverse momentum the first term in 
Eq. ( |45| ) gives an IR finite contribution. The term in the bracket in Eq. ( ^5| ) is also IR finite. However, the two 
delta functions give separately IR singular contributions [from integration over After integrating over s [with 

for the model based on the Borel transformation or 9{so — s) for that suggested by the local duality], these 
become proportional to the valence qq and qqg wave functions, respectively. Combining the IR finite contribution 
from the pure qq sector (first term in the right-hand side of Eq. (|4^) and a similar one for s — > s') gives 



hard 
r vert 



(s,s') = 2(87rCFas)(2V6)^ 



dxdk± dydl^ 



y[li_ - m_] ■ [U + (1 - 2/)q±] 



51 s- 



x{l 



(k_L - xqi 
x{l — x) 



y{l~y)^ {\±-yq±)' + y{l-y) 



(46) 
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Treating similarly two self-energy contributions for the struck quark, and adding them to the vertex corrections 
gives the total IR finite electromagnetic current qq matrix element 



hard / 
Pcm 



where the form factor S is given by 

dydl± 



x{l — x) 



{2V6)S (s 



{k± - xqj_f 
x{l — x) 



(47) 



S{x,'k±,q±) = 1 - ATrCpas 



167r3 





"(2(l-y) + j/2) 




y) + '^^) 








(U - ym 







(48) 



The contribution from this spectral function to F-^ is thus given by 

F.(O') ^ i^^,Tm(Q') = y"[da;dki]*(2;,ki)5(x,ki,qi)*(x,ki -iq^). 



(49) 



Due to the wave function suppression, both k^/a;(l — x) and (k — a;qj_)^/a;(l — x) are smah as compared to q^ in the 
integral in Eq. (^). Thus, the integral for the Sudakov form factor is dominated by 1^, y — + and Ij^ — yqj_, y ^ 0. 
In these two regions, both transverse 1_l and longitudinal y momentum integration lead to log(q5^). To leading 
order in (k^/a;(l — x))/q\ and ((k — a;q^)^/a;(l — x))/q^^, the result is 



5(x,k_L,q±) « 1 7^ log 



An 



An 



log' 



(k_L - XC[A 



(50) 



Alternatively, we may define Sudakov form factor introducing factorization scale, fi to regularize the small 
virtuality of qqg intermediate state, just as we did in defining hard contribution from the one gluon exchange 
diagrams. We thus rewrite the differences of the two (5- functions in Eq. (B3) as 



6 s 



x{l — x) 



1-e 



x{l — x) 



e 



x{l — a;) 



x{l — x) 
S { s 



ki 



(l-x) y{l-y){l-x) 



(51) 



with Q defined in Eq. (p3|), and similarly for the s'-dependent part. The hard contribution which defines perturba- 
tive Sudakov form factor comes, as before, from the first term on the right-hand side. Collecting hard contribution 
from vertex and self energies results in a contribution to F^^ given by 



-J^^i^, k^)5(x, k^, q^; m)*(x, k^ - xq^). 



(52) 



where S{x,'k±,c[±; fi) is now given by Eq. (|4^ ) with k^/x and (k^ — xci±)^/x in the denominator replaced by 
(1 — x)fi^. For ^ 0^ << qi the /^-dependent Sudakov form factor is thus given by 



S[x, ki, qx; ^) w 1 log 



2 



2n "'° \{l-x)^i'^^ 

The full contribution from the diagram on Fig. lb including self energy corrections becomes 



(53) 



-J^■^{x,k^)Six,k^,^r,^im{x,k^-xql_) + J ^g^3 



^g{x, kr,y, l±)T^r''^{x, ki - xq^) + ^{x, k^)r;^-'-^'fg{x, ki - xqr,y, l± - yq 



haidv 



'^g{x,'k±;y,l±)I^^g{x,'k± - a:qi;y, Ij 



(54) 
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Introduction of the factorization scale modifies the Tgi amphtudes according to 



2 ' D, 



, (1 = 1,2) (55) 



and the modification of the nonvalence wave functions '^g — )■ comes from the subtraction of the IR singular 
part of the valence sector, 



%=^g- V^^Al-OiP.)]^^ 

Even though each individual term in Eq. (^) is /i-dependent the entire sum is //-independent. For /i'^ ~ ~ 
(k^/x(l — x)) « the second term in Eq. ( ^6|) strongly reduces the nonvalence amplitudes, and the vertex 
correction becomes approximately given by the valence hard contribution alone, i.e. F^'^^^{Q^,fi ~ f3). Also in 
this case the two expressions for the Sudakov form factor, Eqs. (^) and (^3|) give almost the same result. 
In the ladder approximation Sudakov logarithms exponentiate and S becomes, 

S.e.p(_gc,l„g»^^)^ (57) 

To get rid of the fi dependence in this case one should sum to all orders in a" contributions from qqg^^ wave 
functions whose IR singularities cancel by mixing with nearby Fock sectors. However, as discussed above for the 
n = 1 case, if /i ~ this residual contribution is expected to be small. 

IV. NUMERICAL RESULTS 

The numerical results are summarized in Figs. 2-4. They are all given for fixed (3 = 0.4 GeV. 

In Figs. 2a, b, corresponding to = 0.25 and 0.4 respectively, the triangles represent the contribution of 
the lowest-order diagram Fig. la, i.e., the purely soft contribution determined by the Gaussian wave function of 
Eq. (p5|). At high Q^, the function F^{Q'^) falls off like For comparison, in Fig. 2a we also show the results 

of calculation of the soft contribution within the QCD sum rule method |l^ (dashed line). The set of circles 
represents the sum F^ + F^ + F^ of the lowest-order contribution and the 0{as) correction to the electromagnetic 
vertex F^ ( |29| ) accompanied by the quark self-energy correction F^ (jsj). As shown in Sec. Ill, the exponentiation 
of the leading double-logarithmic terms suppresses vertex correction and it can be approximated by F^^^^ given 



in Eq. (52). Using S from Eq. ( |57| ) with /t^ = 2/3^, we obtain the curve depicted by the squares. Finally, the solid 
line combines F^'^^^ with the remaining contribution from the one gluon exchange diagrams Figs. lc,d, the 
remaining half of the struck quark self-energy diagrams and the spectator's self-energy diagram. 

Fig. 3 gives the comparison of the form factor calculation for two model wave functions discussed in Sec. II. The 
upper and lower dashed lines represent the soft contribution from Eq. ( p^ for the wave functions given by Eq. ( p5| ) 
and Eq. (p^), respectively. The upper and lower solid lines are the respective results for (i.e., full calculation 
to 0{as) for gluon exchange diagrams and Sudakov-exponentiated vertex correction; = 0.25 was used). The 
overall model dependence is seen to be rather small. As discussed in Sec. II, the one gluon exchange diagrams are 
responsible for the enhancement of the hard contribution at high Q^. 

Inspecting the results shown in Fig. 2a (2b) we may conclude that for '-^ few GeV^ the dominant effect 
comes from the soft region. The next-order contribution is ^ 7% (~ 10%) at '--^ 1 GeV^ and increases to about 
20% (30%) at Q2 ~ 4 GeV^ A large fraction of the enhancement fro m one gluon exchange diagrams is however 
IR dominated. This part of gluon exchange should therefore be considered together with the soft wave function 
contribution and separated from the hard scattering amplitude represented by hard gluon exchange. This effect 
is illustrated in Fig. 4. The contributions from hard gluon exchange given by T^^^^ is plotted together with the 
remaining contributions from the one gluon exchange diagrams. The dashed line is the result of the asymptotic 
formula (cf. Eq. (||)) with the asymptotic distribution amplitude given by Eq. (^). The solid line which approaches 
the asymptotic result from below is the contribution from T^^^{^ with /i^ = 2/3^ . The other solid line is the leftover 
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from the one gluon exchange diagrams. The two sets of points shown by circles and triangles below the asymptotic 
curve represent the hard gluon exchange contribution for = 10/3^ and jj,^ = 0.2/3^, respectively. The upper circles 
and triangles give then the remaining one gluon exchange contributions for these two values of the factorization 
scale, respectively. The net contribution from each of the three sets of points (solid lines, circles and triangles) is 
/^-independent. 

As advocated previously /i ~ leads to the fastest saturation of the asymptotic form factor by the pure hard 
gluon exchange i.e. as increases the choice fj, ^ (3 is optimal for reducing the contribution from nonvalence 
sectors. At ~ 20 GeV^ the hard contribution starts to dominate over the IR sensitive part of the one gluon 
exchange. This is consistent with the results found in Ref. |l^] where sensitivity of the form factor to the cut-off 
imposed on the light cone energy denominators was studied. The discussion on the normalization of the space-like 
pion form factor data can also be found in the recent analysis with the optimal renormalization scale and scheme 
dependence p8[ |. 

V. CONCLUSIONS 

Our main focus in this work was to consistently generate the gluon radiative corrections within an approach 
motivated by the light cone quantization formalism and QCD sum rules. Just as in the latter approach, our 
starting objects are the Borel-transformed Green's functions. To make a link with the light-cone quantization, we 
demonstrated that the action of the Borel transform is analogous to using the Gaussian valence soft wave function 
which was frequently used in the past p9[ |. Such a pion wave function can be made to satisfy the standard light 
cone normalization conditions for the pion decay constant, form factor, etc. The main advantage of the Green's 
functions method is that applying the Borel transformation to the two-loop Green's functions, we generate both 
the radiative corrections to the current matrix element and the non-valence qqg components of the pion wave 
function which are absolutely necessary to secure the gauge invariant and infrared-finite results for the total 0{as) 
corrections to the pion form factor. In addition, all of these wave functions are readily applicable to the light cone 
time ordered pcrturbative expansion of pion form factor beyond the leading twist level in QCD. Further works 
involving the four-point Green's functions and applications to the virtual Compton scattering are in progress. 
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FIG. 3. Wave function dependence of QCD corrections to Q^F^(Q^). 
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